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(57) Abstract: The method of reducing contamination of electrical signals recorded in the presence of repeated interference con- 
tamination comprises obtaining an electrical signal recorded in the presence of a contaminating signal, and detecting a timing signal 
that occurs at a fixed time point during the electrical signal relative to the onset of the contaminating signal. The electrical signal is 
digitized, wherein tiie digitizing begins with the tinting signal. A plurality of digitized electrical signals is analyzed, wherein the elec- 
trical signals are synchronized with respect to the timing signal, to obtain an estimated contaminating signal that is subtracted from 
the digitized electrical signal. This method can be used with electrophysiological signals, such as EEG, ECG, EMG and galvanic 
skin response, and for elimination of noise associated with concurrently used methods such as MRT. The method of noise reduction 
is applicable to recordings of other electrical signals, including audio recordings. 
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METHOD AND APPARATUS FOR REDUCING CONTAMINATION OF AN 

ELECTRICAL SIGNAL 

This application claims the benefit of U.S. provisional patent application numbets 
. 5 60/225,389, filed August 15, 2000, and 60/267,337, filed Febmary 7, 2001, the entire 
contents of each of which are incorporated herein by reference. Throvighout this 
appKcation, various publications are referenced. The contents of these references are 
incorporated by reference herein, in order to describe more fully the state of the art. 

TECHNICAL PTFT H OF THE INVENTION 

10 This invention relates generally to methods and apparatus for signal processing and data 

♦ 

collection that are particularly suited for minimizing artifacts and optimizing signal-to- 
noise in simultaneous recording of electtoencephalographic (EEG) and Magnetic 
Resonance Imaging (MEJ) signals, as well as other environments in which electrical 
signals are subject to repeated interference. The methods of the invention can be applied 
15 to other recordings containing repeated electrical interference, including 

eletromyelographic (EMG), electtocardiographic (ECG) or galvanic skin resistance 
(GSR) signals recording during fMRI and audio recordings or transmissions in the 
presence of 60 Hz noise or electrical transients. 

20 BACKGROUND OF THE INVP^NTTOM 

Electroencephalography (EEG)'and functional MRI (fMRI) induce mutual arrifects 
when recorded concurrently. Electroencephalography (EEG) has been a key tool for 
study of the brain for decades. However, despite its multiple clinical and research uses, 
such as in epilepsy (Ebersole, 1997), sleep staging (Rechtschaffen & Kales, 1968) and 
2 5 psychophysiology, litde is yet known about the underiying generators of EEG activity in 
humans. Functional MRI (fMRI) recorded in concert with EEG can provide a method 
for localizing these sources. By using the EEG signal as a reference for fMRI maps, 
concurrent EEG/fMRI opens a new avenue for investigating specific brain function. 



1 



wo 02/13689 



PCT/USOl/25480 



There remains a need for a system for simultaneous recording of EEG and £MRI, which 
can be xised as a tool to localize sources of the EEG. 

Simultaneoiis recording of EEG and £MRI has proven challenging. Time varying 
5 magnetic (B) fields induce an electromotive force (e jn.f.) in a wkt loop perpendicular to 
the B field direction which, by Lenz's Law, is proportional to die cross sectional area of 
the wire loop and to the rate of change of the perpendicular magnetic field (dB /dt). 
When EEG leads are placed inside the MR scanner, the rapidly changing gradient fields 
and the radio-firequency (RP) pulses required for MRI may induce voltages that obscure 

10 the EEG signal (Huang-Hellinger, et aL, 1995; Ives, Warach, Schmitt, Edehnan, & 
Schomer, 1993). The induced e.m.f. yields currents that can cause heating of the 
electrodes and leads and potentially impart bums to the patient (Lemieux, Allen, 
Franconi, Symms, & Fish, 1997). Motion of the leads themselves within the static field of 
the magnet also induces an e.m.£; even pulsatile motion related to heart beat yields 

15 ballistocardiographic artifact in the EEG that can be of roughly the same magnitude as 
the EEG signals themselves .(Ives, Warach, Schmitt, Edelman, & Schomer, 1993; Mud, 
et al., 1998)* Further, introduction of EEG equipment into the scanner potentially can 
disturb tbe homogeneity of the magnetic field and distort the resulting MR images. 

20 Iq addition to the large artifacts in the EEG caused by high firequency gradient and RF 
pulses, the high pass filters of most EEG equipment lead to long signal recovery times 
once the MR acquisition has terminated (Krakow, et al., 1999). One method used to 
overcome these difficulties in studies of epilepsy has been to monitor the EEG in the 
absence of scanning while the patient is ia the magnet and to then trivet functional 

25 scanning manually afi:er identification of inter-ictal spikes in the EEG record (Krakow, et 
al., 1999; Seeck, et aL, 1998; Warach, et al., 1996). Visual evoked potential has been 
studied using intedeaved blocks of EEG and fMRI, where the same stimuli are presented 
in each block (Bonmassar, Anami, Ives, & Belliveau, 1999). In these methods, EEG and 
fMRI are acquired serially, resulting in protocol limitations and problems with data 

3 p analysis. In the triggered method, relevant changes in the EEG can not be seen during 
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functional scanning. Problems also exist witii non-unifotm MR image contrast, given that 
Tl saturation ^ically does not reach equilibrium until 3 to 4 TRs after initiation of the 
scan (depending on TR and effective flip angle). Most often tibis is handled by ignoring 
images acquired in the first 3-4 XRs, but diis then leads to an inherent time delay in the 
5 functional scanning. This could be mitigated to some degree by using schemes that 

correct for the Tl -related intensity differences based on tihe actual TR (DuBois & Cohen, 
2000; Guimaraes, et al., 1998). In the intedeaved method, in addition to the former 
confounds, the EEG and fMRI can not be compared directly. 

10 SUMMARY OF THE INVENTION 

To overcome the limitations in the prior art described above, and to overcome other 
limitations that will become apparent upon reading and understanding the present 
specificatioxi, the invention provides a method of reducing contamination of electrical 
signals recorded in the presence of repeated interference contamination. The method 
15 comprises obtaining an electrical signal, wherein the electrical signal was recorded in the 
presence of a contaminating signal, and detecting a timing signal that occurs at a fixed 
time point during the electrical signal relative to the onset of the contaminating sigoaL 

The method further comprises digitizing the electrical signal, wherein the digiti^^tig 
begins with the riming signal. A plurality of digitized electrical signals is then analyzed, 

2 0 whereia the electrical signals are synchronized with respect to the timing signal, to obtain 
an estimated contaminating signal The estimated contaminating signal is subtracted 
firom the digitized electrical signal, thereby reducing contamination of the electrical 
signal In a preferred embodiment the analysis to obtain the estimate of the 
contaminating signal comprises averaging the electrical signals. In some embodiments, 

25 the analysis to obtain the estimate of the contaminating signal comprises calculating a 
weighted average of the electrical signals. The estimate of the contaminating signal can 
be biased towards recent events, for example, by adding the nth electrical signal to a 
scalar multiple, w, of the prior estimate of the contaminating signal and dividing this first 
sum by a second sum obtained by adding the series l+w^+w^+w^ ,.. + ... w». The 
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estimate of the contaminating signal furthet can be multiplied by a scakx ptiot to the 
subtracting step. 

The method is particulady suitable for electrical recotdings which comprise an 
electtophysiological signal, such as an electtoencephalogiaphic recording, an 
5 electromyelographic recording, an electrocatdiogtaphic recording or a measure of 
galvanic skin resistance. The method is applicable as well to other types of electrical 
recordings, including audio recordings. In some embodiments, the interference 
comprises interference adsing &om inductively coupled magnetic fields. The 
interference can also comprise interference arising from alternating current (AC) line 
10 noise. 

One ad\^ntageous feature of the method of the invention is that the digiti2ing can be 
performed at a sampling rate below the Nyquist rate for the contaminating signal. In one 
embodiment, the electrical signal obtained is passed through a low pass filter prior to the 
digitizing, at a fiequency of approximately one half of the fi:equency at which the 
15 electrical signal is sampled. For example, the low pass filter may pass signal firequencies 
of less than about 200 Hz. 

The method can be performed concurrently with Magnetic Resonance Imag in g of the 
subject In one embodiment, the electrical signal comprises an electrophysiological signal 
and the contaminating signal comprises gradient activity. Examples of a contaminating 
2 0 signal include radio firequency transmitter activity. In a preferred embodiment of the 
method, die digiri7ing is performed at a rate of about 200 to about 5000 samples per 
second. The digitizing can be performed at rates below 200 and above 5000 samples per 
second, with representative rates including 100, 250, 500, 1000, 2000, 3000, 4000 and 
6000 samples per second. 

2 5 The invention additionally provides a metiiod of removing a DC offset from the 

electrical signal by analog subtraction prior to the digitizing. Preferably, the DC ofifeet is 
measured and subtracted from the electrical signal using a difference amplifier. In one 
embodiment, the DC offset is measured by analog to digital conversion, and averaged 



4 



wo 02/13689 



PCTAJSOl/25480 



ovet a time pedod long compated to the lowest ficequendes of interest in the electdcal 
signal An example of such a long time period is approximately 10 times longer t^^^ the 
lowest ftequendes of interest in the electrical signal For example, where lowest 
ficeqiaendes of interest are approximatdy 3 Hz, the time period is about 30 seconds. In 
one embodiment, the analog subtraction comprises converting the averaged signal to an 
analog voltage and electrically subtracting the averaged signal £tom the dectrical signal 
through differential amplification. In another embodiment, the DC offset is measured in 
an analog integrator having a time constant long compared with lowest ficequendes of 
interest in the signal 

The method of the invention is useful for electrophysiological recordings, such as in an 
dectroencephalogram that is recorded concurrentiy with magnetic resonance image 
acquisition. In a preferred embodiment, the dectcophysiological recordings are used to 
inform interpretations of magnetic resonance images. The dectcophysiological 
recordings can be used in a statistical analysis of change in intensity of the magnetic 
resonance signal 

The method can further comprise determining a correlation between change in intensity 
of the magnetic resonance signal and a feature of the dectrophysiological recording. The 
correlation can be used to make statistical images, or image maps, that represent an 
association between the dectrical signals and the intensity of the magnetic resonance 
signal intensity. In one embodiment, die feature of the dectrophysiological recording 
comprises a time course of signal intensity change in defined fi:equency bands contained 
in the dectcophysiological recording. 

The defined frequency bands can be sdected to correspond to standard ranges used for 
clinical interpretations of the dectroencephalogram. Representative standard ranges are 
selected from the group consisting of firom 0 to approximately 4 Hz (the Ddta band), 
fi:om approximatdy 4 to approximatdy 8 Hz (the Theta ban<^, from approximately 8 to 
approximatdy 12 Hz (the Alpha banc^, firom approximatdy 12 to approximatdy 30 Hz 
(the Beta band), and firom approximately 30 Hz and greater (the Gamma band). 
Typically, the frequency bands in this context will not extend beyond 300 Hz. In one 
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emboditneaty tbe method fiitthet comptises convolving die time coutse of tbe 
electrophysiological signal widi an estimate of die magnetic resonance hemodynamic 
impulse response function. In this embodiment, die time course of die 
electrophysiological signal is suitably conditioned to more accurately reflect the 
5 anticipated time course of the magnetic resonance signal change. 

The invention additionally provides a metiiod of reducing magnetic interference during 
electrophysiological recording from a subject by measuring an electrical potential 
difference between a pair of electrodes, wherein the pair of electrodes communicate with 
a differential amplifier via electrical connections, die mediod comprising twistiog the 

0 electrical connections togedier, thereby rcdudng magnetic interference. Also provided is 
a method of reducing magnetic interference during electrophysiological recording from a 
subject by measuring an electrical potential difference between a pair of adjacent 
electcodes, wherein each electrode comprises two leads, the method comprising twisting 
each lead together with a lead of an adjacent electrode, tiiereby reducing magnetic 

5 interference. In one embodiment^ the electrophysiological recording comprises an 
electroencq>halographic recording. The method can be performed concurrently witb 
Magnetic Resonance Itnaging of the subject 

The invention frtrther provides an apparatus for processing digitized electrical signals in 
the presence of a repeated contaminating signal. The apparatus comprises a signal 
0 processor adapted to receive a recording of an electrical signal; a detector adapted to 
detect a timing signal that occurs at a fixed time point during an electrical signal relative 
to the onset of a contaminating signal; a signal accumulator to contain the estimated 
contaminating signal; and a processor adapted to subtract averaged waveforms from an 
electrical signal. The signal accumulator can be, for example, a signal averager. 

5 BRIEF DESCRIPTION OF THE FIGURES 

Figure 1 A. Digital photograph of chained bipolar dual-lead dress. Leads of consecutive 
electrodes are twisted together to reduce scanner arti&ct 
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Figute 16. Schematic diagram showing electrode comiectot^ 2 attached to the head of a 
subject 1 using standard electrode geL Each of die connectors is attached to 2 electrical 
wires (tTpicaUy constructed of carbon fiber material, which reduces txiagnetic atti&cts). 
The wires firom adjacent electrodes are twisted tightly together in paits 3, where each of 
5 the two wires ftom a single electrode is twisted together with a different neighbor. The 
electrode pairs are dien presented to the input of a differential amplifier 4 \diere the 
electrical potential difference is amplified to form die dectroencephalog^iam. The inputs 
to the amplifiers are bridged sudi that the paired leads enclose a complete loop, thereby 
minimiring additional differential potentials between amplifiers. The wires are drawn in 
10 thick and thin lines for clarity only. 

Figure 2A. Schematic diagram showing how dual lead electrodes allow each bipolar pair 
3 to be twisted together for their entire length, sending signal directly to local differential 
amplifiers 4. 

Figure 2B. Schematic diagram showing how twisting of leads leaves only small loops at 
15 the head in which e.m.f. can be induced. Current induced in lead twists by motion and 
gradient switching will be self canceling. 

Figure 3. Diagram of an EEG'data pathway. EEG signal is fed to a local differential 
amplifier 4, digitized and then sent out of the scanner room 7 via optical fiber for real 
time display and off-line analysis. 

2 0 Figure 4. Traces illustrating the ballistocardiogram subtraction algorithm, shown using 
data collected on a normal volunteer. A) A s^ment of the subject's QRS wave was 
correlated to their EKG, and peak correlation values initiated a trigger (shown by the 
vertical dotted lines). B) Tri^er to trigger segments were then averaged, weighting 
segments less and less by temporal displacement fi:om the nth (shaded) segment. Q The 

2 5 averaged data was subtracted firom the raw EEG to yield ballistocardiogram^firee EEG. 
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Figure 5. EEG of a phantom lecorded inside the MR scanner using twisted and 
untwisted leads, showing recordings both in absence of scanning and during EPL With 
or without scanning, EEG recorded using the untwisted leads was significandy noisier. 

Figure 6. Raw EEG of a normal volunteer recorded inside the MR scanner using 
5 untwisted (top) and twisted (bottom) leads. EEG recorded with the untwisted leads was 
significantly noisier. Ballistocardiogram is visible in the twisted lead data, but is reduced 
compared to the imtwisted leads. 

Figure 7. EEG recorded on a volunteer during £MRI, before (above) and after (below) 
post-processing to remove the ballistocardiogram. 

1 0 Figure 8, Power spectrum of EEG recorded simultaneously with £MRI, in steps equal to 
the TR of 2.5 seconds, showing expected increases in the a^ha band (8-12 Hz) when the 
subject's eyes were dosed. 

Figure 9. Graph depicting how errors in gradient noise cancellation will occur when the 
sampling and gradient activity are asynchronous. In this Figure, the sampling used to 
15 create an error estimate (open circles) has drifted by approximately 200 fis compared to 
the current sample (closed circles). Subtracting the error estimate actually increases the 
. residual error (by more than 17%), as indicated by diamonds. 

Figure 10. Depiction of timing for the gradient echo EPI pulse sequence used in Figures 
11-13. 

2 0 Figure 1 1 . Top: Graph representing uncorrected and corrected signals obtained during 
imagbg, using 10 kHz sampling. BoUom Enlarged (lOS) views of the 25 msec periods 
indicated in dashed lines (uncorrected signal in dotted lines). Note that the artifact 
^ suppression varies from cycle to cycle, as a result of phase errors. 

Figure 12. Graph showing simulation of the effects of sampling rate on the efficiency of 
25 gradient artifact suppression. The trace at the top is actual EEG data recorded at 10 kHz 
during an echo-planar imaging sequence. The three traces below it are the difference 
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between the ojdginal signal and a sample la^d at 100 )is, 200 |xs and 400 fxs as indicated 
(the wotst case ettois fot 10 kHz, 5 kHz and 2500 Hz sampling). The graph below 
shows a detail of the pedod dudng which echo-planar readout occurs. 

Figure 13. Graph similar to that shown in Figure 4, depicting a simulation that shows 
5 the effects of timing errors when sampling at 200 Hz - well below the Nyquist feequency 
for the gradients. The waveforms are clearly undersampled and therefore appear quite 
different firom the prior Figure. The magnitudes of the residual artifacts, however, are 
very similar. The three lower curves show the artifacts that remain after correction with 
timing errors of 100, 200 and 400 ^s. 

10 Figure 14. Schematic representation of a low-cost offset nulling differential amplifier 
circuit for use in fMRI. Power supply connections are omitted for clarity. All capacitor 
values are in microfarads. 

Figure 15. MR gradient activity recorded with triggered 200 Hz sampling using a 3s 
repetition time (IR) and nineteen slices. The uncorrected signal is shown at top for a 
15 single TR. In the middle is shown the average of 30 TR periods, and at bottom, the 
difference between the uncorrected and averaged signals ("'corrected"). 

Figure 16. Human EEG data collected during echo-planar functional imaging. The 
uncorrected data appear at bottom. Above them are corrected records from twenty 
successive TR periods. 

2 0 Figure 1 7A. Graph showing energy as a function of frequency and time, derived from 
EEG data acquired during scanning. 

Figure 17B. Graph showing estimated fMRI activation time course for the EEG data 
appearing in Figure 9A. For clarity, only the Alpha (solid line) and Theta (dashed line) 
bands are shown. 

25 Figure 18A-C. Functional MRI statistical maps of signal change correlated with spectral 
energy at each of five frequency bands (18A left, Delta; 18A right, Theta; 18B left. Alpha; 



9 



wo 02/13689 



PCTAJSOl/25480 



18B right. Beta; 18C, Gamma), e3q>tessed as coefficient of cottelation. Note that the 
color scale for the lower ftequerides (18A) is different, as the correlations were overall 
lower. All five images were calculated fiom the same 4:30 (min:sec) acquisition taken 
with the subject at rest with eyes open. 

5 Figure 19. F\inctional block diagram of MRI-compatibleEEG amplifier. 

Figure 20. Alternative differential amplifier circuit for use with £MRI. 

Figure 21 . Flow chart illustrating method of signal correction by artifact reduction. 

Figure 22. Schematic illustration of elements of an apparatus of the invention. 

DKTATTED DR SC^IIPTIO N 

10 The invention is based on the discovery that contamination of a digitelly encoded 
electrical signal can be reduced significantly by making use of a timing signal that is 
associated with the onset of a repeated contamination signal. Such a timing signal can be 
used to align the digitization of repeated contamination signals for detertnining an 
estimate of the contamination which can then be subtracted firom the electdcal signal. 

15 This method is particularly suited for use with electrophysiological signals, such as EEG, 
ECG, EMG and galvanic skin response (GSR), and for elimination of noise associated 
with concurrently used methods such as MRI. Although the examples described in detail 
herein address the appEcation of the method to recording an EEG in the presence of 
fMRI, those skilled in the art of signal processing will appreciate that the method of 

20 noise reduction is applicable to recordings of other electrical signals, including, for 
example, audio recordings, wherein it is desired to reduce or eliminate one or more 
sources of repeated contamination. 

All scientific and technical terms used in this application have meanings commonly used 
25 in the art unless otherwise specified. As used in this application, the following words or 
phrases have the meanings specified 
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As used herein, "electrical signals synchronized \7ith respect to a timing signal" means 
that data corresponding to electrical signals recorded over a period of time are aligned so 
that a timing signal that occurs within each of the electrical signals recorded over time is 
superimposed by the alignment When an average is calculated of the ^ectrical signals 
5 superimposed in this manner, the timing signal, as well as any other signal that recurs at a 
fixed time point relative to die timing signal, will be enhanced relative to non-recurring 
signals. 

As used herein, "twisted" means united by winding, intertwining or coiling. A pair of 
electrical connections between a pair of electrodes and a differential amplifier is 
10 sufficiendy twisted if the enclosed magnetic fields are substantially reduced 

As used herein, "a" or "an" means at least one, unless the context cleariy indicates 
otherwise. 

Methods 

The invention provides a method of reducing contamination of electrical signals 
15 recorded in the presence of repeated interference contamination. The method comprises 
obtaining an electrical signal, wherein the electrical signal was recorded in the presence of 
a contaminating signal, and detecting a timing signal that occurs at a fixed time point 
during the electrical signal relative to the onset of the contamixiating signal The method 
fiirther comprises digitizing the electrical signal, wherein the digitizing begbs with the 
2 0 timing signal. A plurality of digitized electrical signals is then analyzed, wherein the 
electrical signals are synchronized with respect to the timing signal, to obtain an 
estimated contaminating signal The estimated contaminating signal is subtracted firom 
the digitized electrical signal, thereby reducing contamination of the electrical signal 

In a preferred embodiment, the analysis to obtain the estimate of the contaminating 
25 signal comprises averaging the electrical signals- In some embodiments, the analysis to 
obtain the estimate of the contaminating signal comprises calculating a weighted average 
of the electrical signals. The use of weighted averages can serve to achieve an adaptive 
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arti&ct reduction. The estimate of the containinating signal can be biased towards recent 
events, for escample, by adding the ndi electrical signal to a scalar multq>le, w, of the prior 
estimate of the contaminating signal and dividing this by the sum of the series 
l+w^+w^+xi^ ... + ... The estimate of the contaminating signal can be multiplied 
5 by a scalar prior to the subtracting step. 

Figure 19 is a general functional block diagram of representative analog electronics useful 
in the method of the invention. The method performs optimally with an adequately 
linear electronic signal. With reference to Figure 19, signal is applied differentially at the 
input terminals 42. Using matched passive components for each lead 44, contaminating 

10 signal firom sources such as radio firequency is attenuated before being differentially 
amplified using standard components 46. The differential amplifier 46 is commonly 
provided with an offset reference input, such that voltages appearing at this ti*rtm'fiq1 are 
subtracted firom the output Using a sample and hold device 52, the amplified DC offset 
potential, derived fixsm the inputs, is first measured, then applied to the differential input 

15 amplifier 46. To attenuate sources of contamination that are above the highest 

ficequendes of interest in the signal, an active Ipw pass filter 48 is supplied. The filtered 
signal is buffered in the output amplifier 50 before being made available to the digitizer 
circuit A switching means 54 is provided to the sample and hold drcuit to allow 
detection of the input DC offset at any time desired. 

2 0 The circuit modeled in Figure 19 is an example that is especially smted to the problem of 
recording electrophysiological signals during Magnetic Resonance Imaging. The 
functional logic of this diagram is shown in somewhat more detail in Figure 20 which, in 
addition, illustrates a means of transmitting the EEG data out of the MRI suite. One 
skilled in the art will see immediately that many different circuit topologies are possible 

25 that will accomplish essentially the same function for this or other applications. For 

example, in some embodiments it may be desirable to apply a DC offset after the high 
pass filter (as indicated in the schematic of Figure 14), or to avoid this step altogether, if 
significant DC offsets are not present in the signal, Furdier, in cases where there exists 
adequate headroom for the differential amplifier 46 to remain in its linear range for all 
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expected inputs, it may be desirable to substitute a bandpass filter for 48. Hie passive 
filter 44 may be eliminated when any signals that reach the input to the input differential 
amplifier 46 are of suitably low amplitude to avoid saturation effects. 

With reference to the diagram shown in Figure 20, EEG signal from the subject 1 is 
5 carried via twisted pair leads 3 to die input of a battery powered head amplifier. The 
inputs include RF attenuation via a series iaductance and a parallel capacitor. These 
inputs are mixed in die inputs to shield driver 64, whose output is applied to a concentric 
shield 56 surrounding the twisted leads 3 and connected to the subject 1. The inputs are 
coupled to a differential amplifier 58, whose output is applied to an isolation amplifier, 

10 such as the IS0122 firom Burr-Brown Corporation. This device, and a second similar 
device 62 provide electrical isolation to the subject and an added safety factor for EEG 
recording. The output £com isolation amp 60 is sampled by an analog to digital 
convertor 66, whose digital output is stored in a latch 68 and converted to an analog 
voltage by digital to analog convertor 70. This output is applied as a DC correction 63 to 

15 differential amplifier 60 after electrical isolation by isolation amplifier 62. 

The output fcom this head amplifier is presented to a low-pass filter to attenuate signal 
outside of the desired range of the EEG signal. This single channel output may be 
multiplexed with the outputs of other similar amplifiers by analog multiplexor 74 clocked 
by hardware timer 78. The output of the multiplexor 74 may be converted to an optical 
20 signal by optocoupler 76 and transmitted by pptical fiber, together with the dock signal 
to a second optocoupler 82 that is located outside of the MR shielded room through a 
penetration panel 80. De-multiplexor 82 is used to separate the signals firom multiple 
amplifiers and its output may be presented to a differential line driver 88 for transmission 
of the signals over long distances to an analog to digital convertor for later processing. 

2 5 Figure 21 is a flow chart illustrating an application of the method of the invention. The 
raw digitized signal, containing both die desired signal and contaminating artifact are 
shown as (1) in Figure 21. It is strongly preferred, for the invention to work optimally, 
that the signal in (1) be faithfully (linearly) recorded and that the digitization be timed 
with adequate precision to the artifact To produce the corrected signal, an estimate of 
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the ard&ct (2) is simply subtracted firom the taw signal Optionally, the estimate of the 
atti&ct may be multiplied by an ampUtade constant, to account for differences in the 
magnitude of the coupling of the ard&ct to the signal recording system. 

The artifact estimate can be computed as follows: Each time a new raw sample (1) is 
5 available, it is added to the current estimate of the artifact (2) which has been multiplied 
by a scalar amplitude constant, w (4). This summed signal is then divided by the sum of 
the series 1 + w + w^ + w^ + . . . res\alting in a new representation of the artifact (5), 
which then replaces the value used for (2). When w is a number less than 1, the process 
results in a leaky average, where more recent signals have a greater influence on the 
10 estimated arti&ct (2) than do less recent signals. In this way, the system adapts to slow 
changes in the artifact, if necessary. For the purposes of this disclosure, we have called 
this implementation a 'leaky averager." 

More formally, if Rln] is the n'^ raw signal collected, A[n] is the estimated artifact for 
coUectLon and S[n] is the artifact-free signal: 

B[«]+w4«-lj 

J „ 

i=Q 

n 

=Ml 

n 

and 

S[n]«R[/i]-A[w]. 

In the leaky averager discussed above, the influence of temporally distant data decreases 
with time. It would be possible to determine w adaptively by looking at the history of the 
20 estimated artifact. If the artifact is changing quickly, w should be made smaller (reducing 
the influence of older frames on the current correction). Conversely, w should be large if 



14 



wo 02/13689 



PCT/USOl/25480 



the artifiact is very stable. In the limit, if the artifect cannot change over time, A[»] should 
be the simple average of all of the samples (w=l). If the algorithm is to be used in real- 
time, A[n] will be the average of all samples until time ft. If it is used "off-line" A[n] will 
be the average of all samples both prior to and after time n. 

5 It is possible to determine die value of m adaptively as well, approximating it as the 
amplitude that minimizes the correlation between S[n] and A[«-/]. 

The method is particularly suitable for electdcal recordings which compnse an 
electrophysiological signal, such as an dectroencephalographic recording, an 
electromydographic recording, an electrocardiographic recording or a measure of 
10 ' galvanic skin resistance. The method is applicable as well to other types of electrical 
recordings, including audio recordings. In some embodiments, the interference 
comprises interference arising from inductively coupled magnetic fields. The 
interference can also comprise interference arising £t:om other sources, such as 
alternating current (AC) line noise. 

15 It is particularly advantageous that, in accordance witii the invention, the digitizing can 
be performed at a sampling rate below the Nyquist rate for the contaminating signal In 
one embodim,ent, the electrical signal obtained is passed through a low pass filter prior to 
die digitizing, at a firequency of approximately one half of die firequency at which the 
electrical signal is sampled. For example, the low pass filter may pass signal firequendes 

20 of less than about 200 Hz. 

The mediod can be performed concurrendy with Magnetic Resonance Ttnagifig of the 
subject In one embodiment, die electrical signal comprises an electrophysiological signal 
and the contaminating signal comprises gradient activity. Samples of a contaminating 
signal include radio frequency transmitter activity. In a preferred embodiment of die 
2 5 method, the digitizing is performed at a rate of about 200 to about 5000 samples per 

second. The digitizing can be performed at rates below 200 and above 5000 samples per 
second, witii representative rates includiag 100, 250, 500, 1000, 2000, 3000, 4000 and 
6000 samples per second. 
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The inveation additionally provides a mediod of removing a DC offiset from the 
electrical signal by analog subtraction prior to the digitizing. Preferably, the DC o£&et is 
measured and subtracted from die electtical signal using a difference amplifier. In one 
embodiment, die DC offset is measured by analog to digital conversion, and averaged 
5 over a time period long compared to the lowest frequencies of interest in die electrical 
signal An example of such a long time period is approximately 10 times longer dian die 
lowest frequencies of interest in die electtical signal For example, where lowest 
frequencies of interest ate approximately 3 Hz, the rime period is about 30 seconds. In 
one embodiment, the analog subtraction comprises converting the averaged signal to an 
10 analog voltage and electrically subtracting the averaged signal from the electtical signal 
dirough differential amplification. In anodier embodiment, the DC o£kct is measured in 
an analog integrator having a time constant long compared widi lowest frequencies of 
interest in the signal 

The method of the invention is usefiil for electrophysiological recordings, such as an 
decttoencephalogram that is recorded concurrendy with magnetic resonance image 
acquisition. In a preferred embodiment, the electrophysiological recordings ate used to 
inform interpretations of magnetic resonance images. The electrophysiological 
recordings can be used in a statistical analysis of change in intensity of die magnetic 
resonance signal The method can furdier comprise determining a correlation between 
change in intensity of die magnetic resonance signal and a feature of die 
electtophysiological recording. The correlation can be used to make statistical images, or 
image maps, diat represent an association between the electtical signals and die intensity 
of die magnetic resonance signal intensity. In one embodiment, die feature of the 
electrophysiological recording comprises a time course of signal intensity change in 
defiuQed frequency bands contained in the electrophysiological recording. 

The defined frequency bands can be selected to correspond to standard ranges used for 
clinical interpretations of the electroencephalogram. Representative standard ranges are 
selected from the group consisting of from 0 to approximately 4 Hz (the Delta band), 
from approximately 4 to approximately 8 Hz (the Theta banc^, from approximately 8 to 
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approximately 12 Hz (the Alpha band), from approximately 12 to approximately 30 Hz 
(the Beta band), and from approximately 30 Hz and greater (the Gamma band). 
Typically, the frequency bands for this application will not extend beyond 300 Hz. In 
one embodiment, tbe method further comprises convolving the time course of the 
5 electrophysiologLcal signal with an estimate of the magnetic resonance hemodynamic 
impulse response function. In this embodiment, the time course of the 
electrophysiological signal is sxiitably conditioned to more accurately reflect the 
anticipated time course of the magnetic resonance signal change. 

The invention additionally provides a method of reducing magnetic interference during 
10 electrophysiological recording from a subject by measuring an electrical potential 

difference between a pair of decttodes, wherein the pair of electrodes communicate with 
a differential amplifier via electdcal connections, the method comprising twisting the 
dectdcal connections together, thereby reducing magnetic interference. Also provided is 
a method of reducing magnetic interference during dectrophysiolog^cal recording from a 
15 subject by measuring an decttical potential difference between a pair of adjacent 

decttodes, wherein each dectrode comprises two leads, the method comprising twisting 
each lead together with a lead of an adjacent dectrode, thereby reducing magnetic 
interference. In one embodiment, the dectrophysiological recording comprises an 
dectroencephalographic recording. The method can be performed concutrendy with 
20 Magnetic Resonance Imaging of the subject 

A pparatus 

The invention further provides an apparatus for processing digitized electrical signals in 
the presence of a repeated contaminating signal. The apparatus comprises a signal 
processor 90 adapted to lecdve a recording of an dectdcal signal; a detector 92 adapted 
25 to detect a timing signal that occurs at a fixed time point during an decttical signal 
relative to the onset of a contaminating signal; a signal accumulator 94 to contain the 
estimated contaminating signal; and a processor 96 adapted to subtract averaged 
waveforms frona an electrical signal Representative variations of the apparatus are 
described in Figures 3, 14, 19, and 20. 
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The signal processor 90 adapted to receive a recording of an electrical signal may be, for 
example, an electronic circuit (IC) consisting; of a integrated differential amplifier such as 
a an INA114 from Burr Brown Corp., with a DC offeet reference provided through 
measurement by a sample and hold IC such as die LF298 from National Semiconductor 
5 corporation, a low pass active filter and an output buffer all made using standard 

operational amplifier ICs. In this embodiment, the detector can also include a means of 
analog to digital conversion, such as a National Instruments NI 6031E installed in a 
personal computer. 

The detector 92 adapted to detect a timing signal that occurs at a fixed time point during 
10 an electrical signal relative to the onset of a contaminating signal may be, for example, 
an optoisolator IC whose output is conditioned using an IC such an LN555 from 
National Semiconductor corporation to produces a TTL compatible tri^er signal which 
is then presented to an analog to digital converter such as a National Instruments NI 
6031E installed in a personal computer. 

15 The signal accumulator 94 may be implemented, for example, in software in the 'C* 

piograromiog language as a vector of numbers, or in the National Instruments Lab View 
progr amming language as an array of numbers. Either of these may be executed on a 
personal computer. The signal accumulator 94 can be, for example, a signal averager. 
Other methods, in addition to signal averaging, can be used to generate an estimate of 

2 0 the contaminating signal. 

The processor 96 adapted to subtract averaged waveforms from an electrical signal may 
be implemented, for example, as a processing routine in the *C' programming language or 
the National Instruments LabView progcamming language running on a personal 
computer. Those skilled in the art will appreciate variations on the above examples of 
2 5 apparatus elements that will serve the same processing, detecting and accumulating 
functions in accordance with the methods of the invention. 
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Ov^mew of EEG and fMRI in Simultaneous Recording and Mapp in g 

Electtoencephalography (EEG) is established firmly as a means to probe changes in 
electrical signals recorded &om the scalp diat accompany behavioral tasks, and "as a 
marker for clinical, cognitive or neural states. Determination of the three-dimensional 
5 localization of the EEG signal is ambiguous because the relationship between the actual 
position of multiple electrical dipoles and the distribution of electrical potentials detected 
at the scalp has no unique solution. Functional Magnetic Resonance Imaging (fMRI) uses 
modulations in the magnetic resonance signal that depend on variations in blood 
oxygenation to distinguish brain regions whose activity is increased or decreased with 
10 task demands. The following describes a set of solutions to the technical problems in 
simultaneous recording of fMRI and EEG and shows that data from die two methods 
may be combined to create tomographic images indicating brain regions whose activity 
changes as a function of EEG signal intensity in the classically defined spectral bands. 

Significance and Interpretation of the EEG 

15 Study of the electroencephalogram (EEG) is more than a century old (Caton 1875). The 
phenomenon is highly robust namely that electrical potentials exist at the sur&ce of the 
head that are correlated stcong^ with ongoing cerebral activity and fluctuate with sleep 
stages (Rechtschafifen and Kales 1968; Buchsbaum, Mendelson etaL 1982; Benca, 
Obermeyer etaL 1992), emotional state pavidson, SchaflFer etaL 1985; Davidson 1988; 

2 0 Ekman, Davidson et aL 1990; Lambert and Robertson 1999), attention (Klimesch, 
Doppehnayr etaL 1998; Wrobel 2000), therapeutic drug doses (Loo, Teale et aL 1999; 
Alvarez, Lombardi etaL 2000), traits, such as "aggressiveness" (Fishbein, Heming et aL 
1989) and with circulatiog levels of a wide variety of drugs of abuse (Cezayirli, litde et aL 
1975; Maykut 1985; Tokunaga, Takeichi et aL 1989; Abraham and Duffy 1991; Mannelli, 

25 Janiri et aL 1993; Bauer, Gross et aL 1997). Despite the considerable history and attention 
to this measurable phenomenon, the origin of the EEG signal and the locali2ation of its 
sources (presumably cerebral), is still not known. The situation is somewhat more 
favorable with evoked potentials (EP's, in some contexts known as evoked response 
potentials or ERP's) for which the temporally discrete nature, and the motion across the 
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scalp, of the sur&ce potentials cotnhine to give a better indication of the deep dipole ot 
dholes that generate the signals, but it has been difficult to test direcdy the rektionships 
between scalp BP and brain acdvity, especially for sources significandy below the cortical 
sur£a,ce. For example, attempts to locali2e the generators of the brainstem auditory 
5 evoked response by simultaneous recording with depth electrodes or correlation with 
lesions have been conclusive for only a subset of the waveform components present in 
the signal (Starr and Achor 1978; Chiappa, Gladstone Bt oL 1979; Achor and Starr 1980b; 
Achor and Starr 1980a; Goldie, Chiappa et aL 1981; Cohen and Britt 1982; Chiappa and 
Young 1985). 

10 When Caton first reported on the resting and task evoked electrical activity of the brain 
. in animals (Caton 1875), he was able to determine that there was a transcortical electrical 
potential diat changed during periods of functional activity (sensory stimulation). Caton 
noted that, ^^Feebk cmrents of varying i&vction. . were generally present between different 
points on the cortical sur&ce. Some years later, Berger noted that scalp potentials could 

15 be recorded in humans with properties similar to those of Caton's cortical potentials 
(Berger 1929) and he soon reali2ed that diis electroencephalogram varied according to 
die mental state of the subject (Berger 1930). By 1930 Beiger had described what he 
called die alpha rhjrtimi, being relatively high amplitude oscillatioxis in the range between 
8 and 12 Hz that were associated with drowsiness. It is now accepted that alpha activity 

20 is associated witii a relaxed, awake state, usually with eyes closed. 

EEG is now a routine and essential test in clinical neurology. It provides diagnostic 
information that cannot be gathered through any other commonly obtainable means. 
Indeed, its indispensability derives from this lack of other routine clinical tools to assess a 
broad region of cerebral neurophysiology with high temporal resolution, EEG depicts 
2 5 moment-to-moment changes in cerebral cortical function, and thus is valuable in any 
clinical context where such information would guide medical decision making. Such 
situations are not limited to electrophysiologic abnormalities, as other patiiologic 
processes often affect neuronal function and, thus, impact the EEG. These include 
ischemia, metabolic alterations, mass effect, and infection among others (Matkand 1984). 
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Because it tepreseats abnotmality in neuional diectzophysiology, epilepsy is a common 
clinical ptoblem that wamnts the use of EEG (Engel 1984). Despite a central place in 
clinical neutology that has endured over most of the past century, the basis of EEG is 
still understood poorly. The generators of the potentials that sum to create the waves are 
5 xmderstood better (McNamara 1994). However, die interaction between cell populations 
to create electric fields at die scalp has proven difficult to study. 

By performing fMRI during EEG, one may obtain complementary information and 
greater understanding of EEG and the clinical conditions that EEGs may radicate. It 
appears that solid results that characterize die relationship between scalp potentials and 

10 local brain activity will offer great value in guiding die cMnirian to isolating specific 

abnormalities. In fact, there ate already scattered reports diat interictal spike discharges, 
associated strongly widi clinical epilepsy, might be used in combination with fimctional 
MRI (Warach, Ives et aL 1996; Seeck, Lazeyras et oL 1998; Krakow, Woermann €t aL 1999; 
Patel, Blum etaL 1999; Symms, Men etaL 1999; Schomer, Bonmassar etoL 2000) (or 

15 PET (Henry, Sudierling aL 1991)) to identify surgically resectable lesions , and 
tomographic localization via MEG has been suggested as a means to guide such 
resections (Stefaii, Schneider tf/ aL 1990). The use of fMRI during epileptic seizures has 
also been tested with success by Jackson (Jackson, Connelly aL 1994). 

20 EEG and Tmagin g jg Sleep Disotdets and Staging 

Based on unit activity, stimulation studies, and lesion studies in non-human mammals, 
brain regions active in non-REM sleep include the anterior hypothalamus, dorsal bulbar 
reticular formation, and nucleus of the solitaty tract (Jones 2000). Regions likely 
associated widi REM activity, and possibly the generation of wakefulness, include the 

25 posterior hypothalamus, ventral mesencephalic pons, basal forebrain, and pontine 
reticular formation. Such localized brain activity should also be visible in humans 
through functional neuroimaging. Both PET and SPECT have been used to examine 
regional activity changes (via cerebral metabolism and blood flow respectively) with sleep 
stage (determined electroencephalographicaUy). These imaging studies have broadly 

3 0 indicated tiiat areas thought to be involved in die active generation of rapid eye 
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movement (REI^ sleep are active dudng this stage of sleep and in non-xapid eye 
movement (NREM) sleep. 

Maquet and colleagues (Maquet, Dive et al. 1990) noted that rates of cerebral glucose 
5 metabolism (rCMRGlc) during NREM are lower, ovetall, compared to those during 
wakefulness (most notably in the thalamic nuclei), and that REM sleep rGMRGlc is 
comparable to the awake state. Furthermore, the greater the depth of the NREM sleep 

the greater the amount of cortical synchrony), the lower the rCMRGlc (Ingvar, 
Baldy-Moulinier et aL 1965; Madsen, Hohn et al. 1991; Maquet, Dive et aL 1992). 

10 Especially intriguing is the observation that changes in cerebral blood flow show 
substantial regional heterogeneity. REM-associated increases in rCBF have been 
observed in die pontine tegmentum, thalamus, limbic areas, cortical areas (notably the 
anterior dngulate cortex), and visual association areas, with a decrease in the dorsolateral 
prefrontal cortex, parietal cortex, posterior cingulate cortex, and precuneus (Madsen, 

15 Hohn et aL 1991; Madsen, Schmidt et al. 1991; Maquet, Peters et al. 1996; Nofizinger, 
Mintun et al. 1997; Braun, Balkin et aL 1998). Interestingly, rCBF increases have been 
observed in the exttastriate visual areas, though not the primary visiial cortex which may, 
as has been hypothesized by diose authors, be indicative of some sort of visual memory 
activation during sleep (Braun, Balkin et aL 1998; Maquet and PhilEps 1998; Maquet 

20 1999). 

As an example of the relevance of imaging during sleep, cxirrent studies suggest a role for 
sleep in memory consolidation, based on changes in sleep blood flow as a function of 
daytime activities (Maquet, Laureys et aL 2000). There are also important correlations 

25 between sleep and a variety of psychiatric disorders (Benca, Obermeyer et al. 1992) that 
might expose sleep physiology as a marker for these problemis. But, due to the coarse 
temporal resolution of PET and SPECT (Nofzinger, Mintun et aL 1998) and the 
relatively rapid changes in brain activity diaring sleep, neither of these imaging methods is 
suited ideally for this purpose. Very recentiy, there have been reports of the use of £MRI 

30 in assessment of localized signal changes that take place during sleep (Lovblad, Thomas 
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et aL 1999; Home 2000) that indicate an increase in occipital lobe activity and a decrease 
in the frontal lobes during REM, in agreement with the PET findings. Howevei^ such 
studies cannot be considered definitive without the incorporation of sleep staging by 
electroencephalographic means (Lovblad, Thonaas et aL 1999). EEG-£MRI provides an 
5 excellent solution. 

The sleep electroencephalogram can be defined by characteristic firequency patterns, and 
brief electrophysiological phenomena such as k-complexes and sleep spindles. 
Furthermore, during REM sleep, there is a descending suppression of musck tone, 

10 saccadic eye movements, and the loss of the cortical synchrony that is a hallmark of 
NREM sleep. Only with temporal resolving power of fMRI will it be possible to study 
activations associated with these transient events. One application of the invention is to 
assess changes in regional brain activity using fMRI and to correlate such activity to 
classically defined sleep architecture and features. For example, the invention can be used 

15 to seek an understanding the brain activity tiaat undedies the general lack of 

responsiveness to external stimuli, the apparent gating of motor output and dream states. 

Functional MRI (£MRI) is now an established method for the locali2ation of focal areas 
of brain activity, chiefly in humans (Cohen and Bookheimer 1994). Although it is 
assumed that the flNffil signal arises firom local changes in blood oxygen content (Ogawa, 

20 Lee etaL 1990a; Ogawa, Lee etaL 1990b; Kwong, BelUveau etaL 1992; Ogawa, Tank etaL 
1992), this theory has not been subjected to extensive direct testing, and the mechanism 
of coupling between neuroelectncal activity and MRI signal changes is still the subject of 
speculation. Nevertheless, the observed areas of signal increase correlate well with both 
extensive literature on neurophysiology and, more recently, with electro-cordcography 

25 derived of humans in surgical settings (Schulder, Maldjian HaL 1998; Roux, Boulanouar 
etal 1999; Lurito, Lowe ttoL 2000); While there is a growing literature that attempts to 
use the localization power of £MRI to aid in die interpretation of the study of evoked 
responses, there is a striking paucity of reports that attempt to reconcile the findings in 
EEG witii functional MRI. 
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The classical ox postetiot alpha xhythtn is found mostly in ocdpital, panetal,. and 
posteriot tempotal tegioas (Adrian and Matthews 1934), and fitst emerges at about 4 
months of age as a 4 H2 osdflation, present \5ith eyes closed and blocked \dth eyes open. 
The frequency of dais daythm increases witii age, reaching about 8 Hz by age 3, and by 
5 about age 10, it reaches die average adult frequency of 10 H2 (Petersen and Eeg- 

Olofsson 1971). Aldiough recognized since die beginnings of EEG, litde is known about 
die functional significance of diis alpha rhydun; it reflects essentially a state of relaxed 
wakefidness, and can be used as an indirect measure of brain activation, for increased 
alpha band acthdty is thought to correspond to decreased activity in undedying cortex 

10 (Shagass 1972). Thus, decreased a^ha activity, or stimulus-induced alpha blocking, has 
been termed "event related EEG desynchronization" (Pfurtscheller and Aranibar 1977). 
Davidson and colleagues have shown tiiat alpha asymmetry recorded in anterior regions 
correlates widi emotional reactivity, and that these asymmetries appear trait-like in 
waking and in sleep (Petersen and Eeg-Olofeson 1971). Smdies in animals have 

15 suggested the thalamus as a possible generator of the alpha rhythm (Petersen and Eeg- 
Olofsson 1971). Lopes da Silva demonstrated significant thalamocortical coherences in 
dogs between lateral geniculate nucleus and pulvinar and the cortex (Lopes da Silva, 
lierop et aL 1973; Lopes da Silva, Vos tt aL 1980). Recently, in humans, Lindgren and 
colleagues showed an inverse correlation between EEG alpha power and thalamic 

2 0 metabohc rate in normal subjects using PET (Lindgren, Larson et aL 1999). 

Challet^es in Comb ining fMT^ T agd EEG 

Even in die best of citcumstances, EEG signals recorded in die clinical environment are 
rdatively noisy. The effective iaput resistances are large, and die signals are smaL As a 
result, Boltzmann noise limits the ultimate signal to noise ratio. Widi a typical 5 MSl 
2 5 input impedance, die Boltzmann noise wiH be approximately 1 .5 jxV (vans = -yl^kTBR , 
where k is Boltzmann's constant, Tis the temperature, £ is die bandwidth and R is the 
equivalent resistance) even over die limited bandwiddi of 100 Hz or so used in the EEG. 
Since the scalp potentials are typically only a few fxV, die signal to noise ratio (SNR) of 
the EEG seldom exceeds 100:1. Because EEG is often analyzed in die spectral domain 
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over nattow bandwidtihs, the effective SNR for detection of band-limited signal (for 
example, alpha intensity) may be somewhat higher. 

Beyond the thermal noise limit, other factors further reduce the SNR of EEG. Cardio- 
electric (EKG) and myoelectric (EMG) signals add contamination^ Corruption by the 
5 EKG is variable across individuals; though always present, it is usually smaller daan the 
EEG, but at times can become comparable in amplitude. The EMG is typically a 
contaminant for only brief periods of head or facial muscle contraction, such as eye 
blink, grimace, etc. The methods of the invention actually serve to reduce the EKG 
artifact that would be present whether or not the subject is scanned in the magnet All 
10 told, these sources of noise reduce the useable dynamic range of the total EEG signal 
dramaticalLy for virtually all purposes. 

Noise in the MRI Envj rnninpnf. 

The noise environment for EEG becomes radically worse when subjects are placed 
inside an MR imaging system. Almost all of the noise sources, however, are coupled 

15 magnetically to the EEG. Several of these are non-biological, including: amplifier noise 
in the shim and field gradient amplifiers; large time-varying magnetic fields induced by 
the field gradients during scanning, and radio-fiiequency signals generated by the scanner 
for magnetic resonance induction. By Faraday's law, the magnitude of the voltage 
induced by these time varying fields is propordonal to the first time derivative of the flux, 

20 and thxis to die amplitude of the magnetic field, its first time derivatrve, and the area 
enclosed by any conducting loop. More specifically; 

em./ -d^ldt, 

where e,m,f. is the induced dectroinotive force, and (Pis the magnetic flux. In MR 
imagers these sources cannot be reduced in any practical way. The MR imaging gradients 
25 in state-of-the-art imaging instruments slew at an extremely high rate; die field gradients 
on a typical scanner operate at 80 T/ sec and are thus major somces of noise; the newest 
generation of MR instruments, with local head gradients, will slew two to three times 
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faster. The RF p\Jses, though only about 50 tnilliGauss, have fundamental firequendes^ at 
the 3 Tesla operatmg field, of 128 MHz. Slewing at 4000 T/s, they too are large sources 
of noise. 

Physiological noise sources are present also. Even small motions of the subject are 
5 coupled to the EEG as the leads move within the large static magnetic field. Not only 
minor fidgeting, but also the motion of the whole body witt each heartbeat (the 
ballistocardiogram), produces signal in the \iV to mV range. 

DC Offsets and Transient Recovery 

The scalp electrical potentials used in the EEG contain both time varying and static (DQ 
10 components. Often the DC signal is much laj^ than the EEG, but it is seldom of 

interest for clinical diagnostic purposes (note reference to true DC offsets, not to slowly 
varying potentials), as it contains essentially no information. However, it does cause 
trouble for the EEG in several ways. Typically, the DC offset increases the dynamic 
range needed to digitize the EEG signal For example, the EEG signal may be only a few 
15 |iV, while potentials of a few mV may exist between electrodes, or as a result of the 

chemical electrode potential. The signal digitization depth will be reduced by the ratio of 
the EEG to the DC potential For example, assuming a DC offeet of 10 mV and an 
EEG signal of 10 pV, the 4096 different levels representable by a 12 bit analog to digital 
converter (ADC) will be reduced to only 4 levels for the EEG. Cleady this loss is 
2 0 unacceptable, as quantization noise will dominate the signal 

For these reasons, conventional EEG amplifiers are equipped with AC-coupled (high 
pass) inputs, usually a capacitor separating the first stage amplifiers firom the input to the 
ADC. The inputs will usually have time constants of several seconds, allowing 
frequencies of 1 Hz or so to pass unattenuated. One consequence of this AC coupling is 
2 5 that it creates a time constant for signal recovery if the input saturates. Because these 
filters must pass very low fcequencies, the recovery time for ibe analog signal to come 
back to the center of its nominal range can be quite long. 
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In £MRI, the recovery time associated with AC coupling is a substantial problem, as the 
gtadient-induced artifacts (tens of mV) can be large enough to saturate the input stages, 
pinning them to the positive or negative supply rails for several msec. When the 
gradients cease, the setding time of the hi^ pass filter greatly oudasts the gradient event 
5 A recent paper studying EEG-£MRI combinations, reported that, '*The EEG could not 
be interpreted during the artifect caused by the excitation pulse, but the recording 
becomes readable in less than 1 second (approximately 100 msec) after completion of the 
BURST" (Hennig and Hodapp 1993; Lovblad, Thomas et al. 1999). This problem can be 
mitigated by tising the methods of the invention, including an artifect-redudng electrode 
10 configuration and an input amplifier (Grass-Tele&ctor) with enough headroom to stay 
out of saturation prior to the high pass stage. As disclosed herein, the invention further 
provides a more complete solution that avoids the hig^ pass filter completely. 

Gradient Noise 

The magnetically-induced gradient noise is of very large amplitude (tnilliVolts in a typical 
15 scanner) as compared to the EEG, especially in the context of echo-planar imaging. One 
group has implemented a correction scheme for gradient artifacts that is siirtilar to a 
scheme that group developed for ballistocardiogram removal (Allen, Josephs et aL 2000). 
Because the fundamental fi:equencies of the gradient activity are much higher, they 
developed special recording hardware that allowed them to use a much higher digital 
2 0 sampling rate of 5 kHz, which they selected as being rapid compared to the nominal 
Nyquist frequency for the gradient waveshapes. Unfortunately this is not sujBScient, as 
sampling at the Nyqxiist rate guarantees only against aliasing of the higher firequencies 
into the pass band, but does not effectively remove the artifact 

For example, assume that there is an undesired signal consisting of a sinusoid at 1000 Hz 
25 contaminating the EEG, which is sampled digitally at 5 kHz. Because the gradient and 
digitization activity are clocked independendy, and are asynchronous, the phase at which 
the artifact is sampled can vary by as much as 2%/S (12^ or 200 pjs at this frequency). 
Over extended sampling penods (typically five minutes or more in an imaging 
experiment), it is likely that the relative timing of the scanner and sampling clocks will 
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differ to tihis degree. Figure 9 shows diis effect In this simulation, the artifact is assumed 
to have been sampled at the points indicated by closed circles. At a later time, the 
sampling has drifted with respect to die gradients by 200 \xs (dashed line, open circles). 
The difference signal is the residual artifact, in diis case just over 17% greater th^ n the 
5 uncorrected signal. To mitigate this problem, the Allen group has adopted a sophisticated 
interpolation scheme that is successful in minimiTing the residual contamination. The 
present invention provides a much more effective approach based on an alternative 
formulation of the digital sampling problem. 

The methods disclosed herein have made it possible to effectively eliminate 
10 contamination of the EEG signal by the most severe sources of noise present during 
MRI scanning in general, and in functional MRI in particular. In one embodiment 
described herein, the method has been used in the construction of tomographic maps of 
brain activity corresponding to the eneigies in specttaUy-defined components of the 
EEG. 

15 Theory 

Digitization 

The artifact ftrom gradient activity is laige and contains substantial energy at hig^ 
fcequendes. Figure 10 shows the timing of a typical echo-planar imaging sequence, as 
used in typical functional studies. The Unes for Select, Phase and Readout indicate the 
20 amplitudes of die three ordiogonal magnetic field gradients used for imaging. The fourth 
line indicates the timing of the radio fcequency chaimel (only the amplitude envelope is 
shown for die RF, as the carrier frequency of 128 MH2 is not visible at this resolution.) 
Immediately apparent is the very large hig^ frequency (1400 Hz) oscillation of the 
readout gradient (shown at half the vertical magnification of the other gradients). 

25 Figure 11 shows raw signal, recorded fi:om an EEG system at a sampling rate of 10 kHz, 
following analog low-pass filtering at 100 Hz. The insets on the bottom of this figure 
show expanded representations (10!X) of the indicated regions of the signals. When EEG 
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data axe acquired dvudng scans, iihe MRI field gradients induce voltages much higher than 
the cortical signal Compatiog this and the previous figure reveals that the magnitude of 
the high fi:equency components is reduced dramatically by the low pass filter, and that 
the residual arti£acts outlast the gradients themselves. 

5 Though the low pass filter provides at least 100-fold reduction in the 1400 Hz 

oscillations, it does not remove the transients as these gradients turn on and off. These 
contain energy at very low ftequencies, as well. The AC coupled input stage, in this 
traditional amplifier, is responsible for the extended ring-down of the artifacts (the 
saturation recovery alluded to previoxisly), although it is much better than the 0.1 to 1 s 
1 0 ring down reported by others (Lovblad, Thomas et aL 1 999), presumably due in large part 
to the attenuation provided by the differential recording apparatus, which helps to 
prevent the amplifiers firom going into saturation. 

When applying cyclic averaging techniques to this data set (Allen, Josephs etaL 2000) as 
shown in Figure 11, they are reasonably effective in attenuating the effects of the low 
15 firequency components, they are largely //^effective at removing the high ftequency 

content. This, as described above, results ficom the asynchronous sampling. The residual 
(worst case) error firom sampling too slowly can be predicted firom the maximxim phase 
shift, tp 

€ = cos(275/? + ^) - C0S(275^) 
= cos(275^)(cos^-l)-sin(275/?)sinp 

2 0 The maximvim phase shift, tp, diat can occur at a given sampling rate is equal to 27cfc/^, 
where is the sampling firequency, and ft is the firequency of the EPI readout 
Comparing tbe residual artifact during the two scan periods (tsvo expanded firames at 
bottom of Figure 11) reveals that the cancellation efficiency is unstable as a consequence 
of the asynchronous timing of the gradient activity and sampling device, which causes 

2 5 tiie sampling offset, cp, to drift over time. 
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As the sampling rate is increased, the cancellation will become mote accurate. Using the 
approximation that for small a, sina ^ a and cosa ^ 1, one can see diat in diis regime sis 
approximately proportional to cp. It follows that if the artifact must be suppressed by 100 
fold, the signal must be digitized at approximately 100*2*7c times the highest ficequency . 
5 of interest, in diis case (the 1400 Hz readout) about 880 kHz/channel, which is 

impractical for reasons of bodi cost and overall data handling. In any case the acceptable 
sampling error is predicted readily by this formula given the low pass filter characteristics 
and the desired final signal to noise. 

The data above were sampled at 10 kHz, well above the 2800 Hz Nyquist criterion for 
10 the high fi:equency components of die signal Figure 12 shows die efficiency of artifact 
subtraction based simply on rapid sampliag. On the top is shown the raw artifact Below 
it is shown the residual artifact that remains after subtraction if the timing of the 
sampling and the scanner have drifted firom synchrony by 100, 200 and 400 psec, 
corresponding to the worst case errors for sampling at 10 kHz, 5 kHz and the 
15 approximate Nyquist rate of 2500 Hz, respectively. It is immediately apparent that the 
residual artifact, after subtraction, is large even with the smallest timing offset. In the 
graph at the bottom of Figure 12, which shows in greater detail only die echo-planar . 
readout segment of that data, one can see that the simple subtraction actually increases 
die magnitude of the artifact, as predicted in the equation above (and in Figure 9). 

20 The efficiency of the subtraction of the gradient artifacts is effectively independent of 

sampling rate, and a repeated single sample, properly timed, can be used to correct for the 
artifact at that time point completely. This general finding is shown in Figure 13 which 
shows the same effects of timing shifts at a low sampling rate of 200 Hz; well below the 
Nyquist rate (the gradient activity is aliased into the digitized EEG signal). As in the prior 

2 5 example, tbe magnitude of the residual artifiicts increases as the sampling is delayed with 
respect to the scan timing. Its magnitude is no worse than that seen with more rapid 
sampling. 

It is clear that if the sampling is timed precisely to the gradient activity, the residual errors 
wiU be elimina ted much more effectively. Perhaps less intuitive is the fact that this does 
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not require Nyquist rate sampling for the artifact feeqaendes. Time (phase) shift and 
firequency can be seen as "duals": Perfect sampling needed for artifect subtraction could 
be achieved by sampling at infinitdy high frequency with arbitrary timing or vnih precise 
synchronization at an arbitrarily low sampling rate. To yield die same 100:1 suppression 
5 that requires 880 kHz sampling, the sample timing would require an accuracy of 1/880 
kHz, or 1.14 fxsec. Ideally, the residual scanner arti&ct should be small compared to the 
thermal noise of the BEG signal With analog filtering and proper recording technique 
(outlined in the examples bdow), die scanner arti&ct is about ten times the amplitude of 
the EEG. Assuming an EEC signal to noise ratio of 100:1, a thousand-fold suppression 
10 is needed in the digital processing, achiercrable with either 8.8 MHz/channel sampling 

(which could come only at tremendous e^ense in the digitization hardware), or with an 
easily achieved 1 14 nsec timing accuracy 

The following examples are presented to illustrate the present invention and to assist one 
15 of ordinary skill in making and using the same. The examples are not intended in any 
way to otherwise limit the scope of the invention. 

Example 1: Argniring - simultaneous EEG and functional MRI 
Methods 

EEG Device and Lead Placement 

20 The EEG device incorporates numerous hardware modifications to reduce artifact in 
concurrent EEG/fMRI, and was provided by Telefactor Corporation (W. 
Conshohocken, PA). Signal is detected from the scalp using silver chloride plated plastic 
cup electrodes connected to a compact magnet-compatible local amplifier (headbox) via 
10 foot carbon fiber leads with a resistance of 1 kQ/ foot This design minitnized both 

2 5 artifact in the MR images and the induction of BP current loops in the lead wires. 

A lead configuration was devised that minimized unwanted current induction by 
recording EEG in a hard-wired montage ming special dual lead electrodes. The lead 
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wires £com each electrode pair were twisted together over their entire leiigth, fortoing a 
chained bipolar montage for each hemisphere (j^2-fS,JS-Pt, t4'i6, i6-o2, o2'p4,p4'c4, d-ff-, 
f^'S^2;jpUp,p43, t345, t5^1, oUp3,p3^3, c3-j3,f3^^f). The dual leads allowed each 
differential pair to be twisted togetiier (Figure lA-B and Figure 2A). The resulting 
5 configuration leaves only small loops at the head in which current can be induced. As 
shown in Figure 2B, currents induced by motion and gradient switching should be self 
canceling, for the current induced in consecutive twists will flow in opposing directions 
in each lead wire. 

The magnet-compatible headbox contains 32 separate channel inputs, each with a 
10 differential amplifier coupled to an RC filter having a time constant of 0.25 msec. 

Sncteen channels were used to record EEG, and two additional channels were used for 
electrocardiogram (EKG). EKG was acquired using a pair of twisted single lead 
electrodes placed above and below the heart on the subject's back. This placement 
minimized lead motion, and thus electrical arti&ct, due to breathing. A scan tri^er 
15 channel was also used to receive a pulse firom the scanner every TR to aid in post- 
processing of the data. The sigoal was filtered in all channels witli a band pass of 0.5-70 
Hz to further attenuate high firequency noise. A single lead was connected at c^^ to 
headbox ground as an added patient safely measure, but was not used as a montage 
reference. 

20 The signal was fed to a battery powered isocoder containing an A/D converter where it 
was sampled at 200 Hz., and the digitized signal was carried out of the shielded magnet 
room via optical fiber to maintain the scanner's electromagnetic isolation. After 
translation to TTL, the data were routed to a Tdefactor Digital BEG (D/EEG) (486 
computer). The EEG data could then be viewed in real time and sent off-line via a 

2 5 lObaseT Ethernet connection to a post-processing and viewing station for further artifact 
attenuation. 
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ScanPrptQCQlt 

All subjects signed a consent form approved by the UCLA Hnman Subject Protection 
Committee prior to MR scanning, which was performed in a General Electric 
(Waukesha, WI) 3T scanner modified for Echo Planar Imaging (EPI) by Advanced NMR 
5 Systems. No visual or auditory stimulation was provided to tihe subjects during 
functional scanning. 

Scout scans of the entire brain were first acquired to localize slice planes parallel to the 
AC-PC line through the occipital cortex. To acquire EEG during functional scanning, 
the scan protocol was then specified to allow windows of readable EEG between 

10 gradient bursts. An EPI sequence was used with TR = 4000 ms, echo time (TE) = 45 
ms, 64 X 64 matrix, 20 cm x 20 cm field of view (FOV), 4 tnm slice thickness, and 1 tntn 
gap to collect 6 slices spaced evenly over the TR period, leaving a 580 ms window of 
readable EEG between each 90 ms period of gradient induced noise (an 87% duty cycle). 
An EPI sequence (TR = 6000, TE = 54, 128 x 128 matrix, 20 cm x 20 cm FOV) was 

15 then acquired coplanar with the fiinctional scans for use as an anatomical reference. 

Artifact Reduction Post Processii^ 

Post-processing and viewing was performed on a Dell Inspiron 3000 Pentium PC. After 
importing data fi:om the D/EEG over the wire, the EEG data was viewed with 
Telefactor Twin sofiware, and processed further using home-built software described 
2 0 below to remove remaining artifacts. The residual artifEict included both noise firom 

e.m.f. induced by the magnetic field gradients (which appeared in the EEG when a slice 
is acquirecQ and ballistocardiogram. The latter occurred in a fairly regular pattern just 
lagging the Q wave of the EKG, but its morphology and amplitude differed in each 
EEG channel 

25 To suppress gradient noise and ballistocardiogram in the EEG record, the gradient noise 
was first removed by blanking the EEG and EKG data for die duration of each MR slice 
acquisition. Following a trigger every TR, the 90 ms data segments containing scan 
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artifiact were teplaced \vith 2eros. Therefore, large deflecttons caused by the MR 
gradients did not corrupt the EEG data during the averaging and subtraction in further 
arti&ct removal 

Next, using a method similar to that of Men et al. (Allen, Polizad, Krakow, Fish, & 
5 Lemieux, 1998), sections of EEG were averaged together following a cardiac trigger to 
yield the ballistocardiographic artifact To identify the initiation of each cardiac pulse, a 
single arti£act-fi:ee QRS wave segment of the subject's EKG, recorded inside the scanner 
when no scanning was taking place, was used as a reference (see Figure 4A). This 
reference wave segment was compared to portions of the EKG data of the same number 
10 of data points as the reference segment that were shifted by one point at a time, and 
calculated a correlation coefficient (CC) for each data pordon. When the CC exceeded 
an empirically selected value (typically 0.7), the peak CC foUowmg this threshold crossing 
was identified to trigger the initiadon of ballistocardiogram averaging and subtraction. 

Figure 4B illustrates the averaging and subtraction algorithm performed on the data in 
15 each EEG channel. Every trigger-to-trigger section of raw EEG data (An) was averaged 
with all preceding sections. Because EEG and EKG should be imcorrelated, this 
method averaged out the EEG signal and left only ballistocardiogram (Bn). Data 
sections were weighted inversely with their temporal displacement from the current 
sample to compensate for slow changes in the ballistocardiographic artifact, calculating 
2 0 the ballistocardiogram in each td^er to tdgger secdon using the weighted average 




with a weighting factor w = 0.9. Thus, the earlier cycles formed an exponentially 
decreasing contribution with a time constant of roughly 10 sections. The averaged wave, 
Bn, calculated separately for each channel, was subtracted from that channel's raw EEG, 
2 5 An, to yield artifact free corrected EEG. 
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Vatiatioxis in section length due to changes in heatt rate were accounted for by averaging 
sections point-by-point The first data point in each section was averaged with the first 
data point in previous sections, the second with die second, etc. Data points at the end 
of longer sections were averaged with corresponding points in other long sections^ and 
5 the weighting fiictors of each point were adjusted accordingly. To avoid subtracting data 
with too few points averaged, die subtraction was performed only if three or more 
points were used in the calcukdon; otherwise, the raw data remained in the final record. 

Q)ianirffti"t» .ation of Noise Reduction 

Noise Reduction due to Lead Dress: Phantom Study 

10 To characterize noise reduction due to the twisted dual-lead dress, scanning experiments 
were performed using a biological phantom - a 9 pound head-siied grocery store 
roasting chicken - and compared twisted vs. untwisted lead arrangements. The leads 
were placed at distances correspondiag approximately to standard international 10-20 
positions on the phantom. The eig^t electtode chained twisted montage detailed above 

15 was placed on the left hemisphere, and on the dgjht was placed a matching dual lead non- 
twisted montage. A scout scan was acquired to position the functional slices to cover the 
phantom. Functional EPI scans were then performed as described previously. 

The EEG data were analyzed to quantify noise reduction due to the lead dress, 
. calculating the loss, in dB, using 30 second data segments as 



20 



dB = 201o g ' 



■2 

twist 



no- twist 



where averages were taken over all twisted and all untwisted channels. Gradient noise 
reduction was calculated by subtracting the square root of the sum of the squared 
voltages without scanniag fi:om that during scatming, and then calculated attenuation as 
above. 
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Noise Reduction due to Lead PfPAS! T^nmfln Study 

To characterize twisted lead noise reduction furdier, the above study was repeated on a 
25-year-old normal male volunteer. Again, the chained twisted lead set was placed on 
the subject's left hemisphere, and the untwisted set on the subject's right The subject's 
5 EEG was recorded inside the MR scanner, both with and widiout scanning. 

Spectral Analysis 

By restdcdng the tuning of the EPI acquisitions to fall outside of the frequency band of 
interest for EEG (e.g., to study alpha activity, die scanning rate must be less than four 
images/second), it was possible to retain useful EEG spectral infotmadon. To illustrate 
10 this, die 25-year-old nortnal volunteer was scanned during diree different tasks known to 
moderate alpha power — a basic eyes open/eyes closed task, a math task, and a 
visualization task. AH studies were performed during functional MRI as described above. 
A baseline eyes open scan was acquired, and then the three tasks were run as follows. 

In the eyes closed task, the subject was given verbal cues during the scan to keep his eyes 
15 open for the first two minutes, closed the next two minutes, then open for the last 

minute. For the next two tasks, the subject was told to keep his eyes closed. Fits^ the 
subject was instructed to count backwards by sevens from the four digit number given to 
him two minutes into the scan. To verify that the task was performed, the subject was 
asked for the number he ended on after the three minutes of counting. Second, the 
2 0 subject was instmcted to visualiae eating his &vorite meal, again two minutes into the 
scan. 

After post-processing the EEG to remove artifact, the alpha power in each TR was 
calculated using software developed in house. With this software, the EEG power in 
user-defined bands with each TR was found using the Fast Fourier Transform. Spectcal 
2 5 power in the alpha band was used as a reference fimcdon to calculate fMRI signal maps. 
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Noise Reductioii due to Lead Dtess: Phantom Study 

EEG tecorded on the phantom using the untwisted leads was substantially noisier th^r y 
diat recotded using their twisted counterparts. Figure 5 shows the recorded EEG in both 
5 twisted and xmtwisted lead channels. When not scanning, the twisted leads reduced 
random noise power by an average of 5.4 dB. During scanning, gradients caused large 
artifacts in both the twisted and untwisted lead sets, but this noise power was reduced by 
an average of 6.3 dB in the twisted leads. 

Twisted vs. Untwisted Leads on a Volunteer 

10 Figure 6 shows EEG data recorded on a normal volimteer inside the scanner when no 
scanning was taking place. Twisting the leads (shown here on the left hemisphere vs. the 
imtwisted on the right hemisphere) reduced noise by an average of 7.5 dB across all 
channels. 

Noise Reducdott Due to Artifact Post-Processing 

15 Figure 7 shows EEG recorded on a vohinteer during functional MRI before and after 
post-processing. The post-processing removed significant gradient and RF artifact, as 
well as ballistocardiogram. 

Spectral Pata 

The data shown in Figure. 8 was recorded on a normal vokmteer during an eyes open, 
20 eyes closed paradigm using simultaneous EEG/fMRI. Here, four slices were acquired 
with a 2.5 second TR. When the subject's eyes are closed, power in the alpha band, 
between 8 and 12 Hz, increases significandy. And as is expected, this alpha signal is not 
present when the subject's eyes are open. 
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Discussion 

Using a combinatioa of analog pre-preprocessing and digital post-processing, this 
method allows one to record clean EEG during functional MR scanning. While the 
EEG is obscured during gradient bursts, it recovers quickly after each slice acquisition. 
5 With this method, then, a trade off must be made between brain coverage in functional 
scanning and the fraction of useable EEG in the data record Consideration must also 
be given to the timing of slice acquisition so that it does not overiap with desired spectral 
frequencies. 

The invention provides a potentially powerful tool for localiadng the sources of vatiovis 
10 EEG waveforms. Because the EEG is acquired simultaneously with fMRI instead of 
serially, it can be used as a direct source for the fMRI reference fimction. In diis way, 
activation maps could be made of any relevant changes in the EEG: spike and slow wave 
patterns in epilepsy, spectral changes, or even event related potentials. 

Example 2: Method for removal of artifacts in simul taneous EEG and fMRI 

15 The strategy outlined below is applicable to any method involving digital subtractive 
noise cancellation, including EEG, EEG/fMRI, as well as any environment in which 
subtracted averages of artifacts are used to reduce noise. 

Methods 

EEG recp rding an d electrode placement 

2 0 The basic approach of differential recording has been described previously (Goldman, 
Cohen et aL 2000; Goldman, Stem et aL 2000). Paired silver electrodes are placed on the 
scalp surface attached with conductive electrode gel and check to ensure that the 
nominal impedance is less than 5 kQ for each electrode. The electrodes are themselves 
attached to two carbon fiber conductors having a distributed resistance of approximately 

25 3 kfl/m. The wires are dressed in pairs witii the leads from adjacent electrodes twisted 
together tightiy to minimize electromagnetic pickiqp. This configuration offers about 6 
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dB attenuation of ard&cts ftom gradient and ballistocatdiogtaphic noise soiitces 
(Goldman, Stem etaL 2000). 

A mplifi rflrint^ 

Fmther imptovements in the overall petf ormance of in-magnet EEG can be achieved 
5 through the use of better analog electronics. Although the ballistocardiogratn contains 
substantial energy in the range of most interest to clinical EEG (from 1 to 50 Hz), the 
other main sources - gradient noise and RF transmit noise — have much higher 
fundamental frequencies. In a typical scanner, for example, the overwhelming majority of 
the gradient-related noise is at a fixed frequency of 1400 Hz, with significant energy 
10 distributions down to 100 H2 or less. The radio frequency energy, of coutse, is at 
radically higher frequencies well outside of the interesting paiss band for EEG. 

Disclosed herein is a very simple circuit that corrects for much of the analog portions of 
the artifact (Figure 14). An initial gain stage features a single pole filter to attenuate the 
large RF signals prior to sampling and provides enough gain to bring the EEG signal into 
15 the mV range without the artifacts causiag saturation. The next stage offers 30 

dB/ octave attenuation at a comer frequency of 200 H2 (so that gamma range EEG is 
readily passed.) 

The AC coupling problem is handled in the fina^ amplification stage, which is arranged 
to include a resettable offset-nulling circuit that stores any DC offset across a low-leakage 

2 0 capacitor. The offset nulling switch can be, e.g., a mechanical switch, or, so that the 

nulling can be performed under digital control as needed, presumably when the software 
detects that die signal is close to digital saturation, a CMOS switch. One can reduce this 
design to use witii a pc and assemble a single channel of this system. Under lab test 
conditions, the circuit is able to hold the DC ofEset to witiiin 1% for 10 minutes at a 

2 5 time, easily meeting the system requirements for MR scanning. 

The AC coupling problem is handled in the final amplification stage, which is arranged 
to include a resettable offset-nulling circuit that stores any DC offset across a low-leakage 
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capadtot. The o£Gset nulling switch can be, e.g., a mechanical switch, ot, so that the 
nulling can be perf otmed under digital control as needed, presumably when the software 
detects that the signal is close to digital saturation, a CMOS switch. This design has been 
reduced to printed drcuitty and a sin^e chaxmel of this system has been assembled. One 
5 can reduce this design to use with a personal computer and assemble a single channel of 
this system. Under lab test conditions, the circuit is able to hold the DC offset to within 
1% for 10 minutes at a time, easily meeting the system reqiiirements for MR scanning. 

Imaging 

Presented herein are examples of activafion mapping data ficom two subjects (one each 
10 for the raw EEG data and for the EEG energy maps). Both subjects were without 

neurological or radiological abnormalities as assessed by a brief neurological inventory 
based on a form developed by the National Institutes of Health NIH and a neurological 
inventory performed by a board-certified neurologist. Beyond lying in the magnet during 
scanning with eyes closed, the subjects performed no explicit cognitive task. 

15 All scanning was performed on a General Electric (Waukesha, Wl^ 3.0 Tesla Signa® 
scanner modified by Advanced NMR Systems (Wilmington, MA) for high performance 
echo-planar imaging (Brady, Cohen et al. 1991; Cohen, Kelley et aL 1996). For testing of 
artifiact rejection, a 19 slice echo-planar data set was collected with a TR of 3 seconds, 
TE =45 ms, 4 mm slice thickness and 3.125 tritn in-plane resolution (64 x 64 scan matri^r 

20 and 20 cm FOV), to achieve appropriate weighting for ''BOLD" contrast effects 

(Ogawa, Lee et al 1990a). For the mapping data, after image-based shimming (Reese, 
Davis et al 1995) and collection of a scout scan, itnaging was performed using a gradient 
echo EPI scan as described above. These data, however, were acquired with a longer TR 
of 4 seconds and only four slice planes. The pulse sequence was modified to include a 5 

25 ^is trigger pulse at the beginning of each TR period, the leading edge of which was used 
for synchroni2ation of the EEG sampling acquisitions. 
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Digitizarion of the EEG 

EEGdata are acquired using a PCI-1200 (National Instruments, Houston, TX) on a pc- 
compatible roicrocomputer. Using LabView (National Instruments, Houston, TX), 
sampling software was developed tiiat responds to tiie leading edge of the scanner trigger 
5 by acquiring a fixed number of samples at a user-specified rate. Specifically, with a TR of 
3 seconds, one can acquire 599 samples at a rate of 200 Hz following each trigger. These 
data are immediately flushed to a file in a background process; to accommodate this 
process one sample point is dropped with each TR. 

Averaging and Artifact Removal 

10 A fcame of EEG data is collected firom each channel with each scan TR, and the fiames 
are averaged separately for each channel to create an accurate representation of the 
scanner artifact One can then remove the atti£act by simply subtracting the averaged 
signal firom its respective channel The next step is to manually inspect the EEG data for 
characteristic artifacts, such as eye blink and facial muscle movement that are recognized 

15 easily by their morphology within and across channels. 

BalHstorardiogram Suppression 

Details of the ballistocardiogram suppression procedvires have been published previously 
by Goldman et al., (Goldman, Stem et aL 2000) and will only be summarized here briefly. 
Because the cardiac-induced motion is neatly repeated with each heartbeat, the resulting 

20 artifact is essentially the same and is superimposed onto the desired EEG. To remove 
this artifact one can detect the heartbeat using the electrocardiographic signal and 
calculate the average balHstocardiogtam over many heartbeats. One can then subtract this 
average firom the EEG signal. Of course, the morphology of the ballistocardiogram 
differs in each of the EEG chatmels as the motion of each lead differs slightiy. Therefore 

25 this averaging and subtraction process is preferably performed separately for each of the 
EEG lead paics. The method is conceptually like that described by Allen et aL, (Allen, 
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Polizzi et aL 1998) but diffets in tbe details by which slow variations in die 
balUstocacdiogtatn ate accommodated. 

Image Processing 

Using software developed in-house, standard FFT methods were used ^ress, Vettetling 
5 et aL 1992) to determine the power separately in five different spectral bands in the EEG 
signal for each TR period Following convention, dato were combined from 0.5 to 4 Hz 
("delta"), 4 to 8 Hz CThela"), 8-12 H2 ("Alpha"), 12 to 30 H2 C*Beta") and from 30 to 
70 Hz ("Gamma"). Using tiiis model of spectral energy as a function of time, one can 
then estimate a prediction of the BOLD signal changes by convolution of the spectral 
10 data with an a priori model of the brain hemodynamic response ftmction (Cohen 1997). 
The convolution introduces a lag in the time course that is thought to represent a 
reasonable estimate of the hemodynamic latencies and, in addition, acts as a low pass 
filter, which tends to reduce somewhat any noise in the EEG data (see Figure 17). 

One can then use scanSTAT (available via die internet at tie URL for the UCLA Brain 
15 Mapping Center, http: / /www-btfliritTiflpping i^rg} to first spatially filter the images and 
then to form statistical maps indicating die correlation between the EEG spectral power 
in each band and the local fluctuations in MR signal intensity. The convention is to 
indicate regions of increasing positive correlation in colors firom red to yellow and areas 
of increasing negative (anti) correlation in colors firom blue to cyan. In particular, regions 
20 of high negative correlation are interpreted as indicating decreased blood flow and 
metabolic activity. 

Results 

Gradient Artifact Suppression 

Tested was the principle method of triggered sampling by recording the analog 
25 waveform tiiat drives the gradient coils, to determine the efficacy of the method 

independent of physiological signal fluctuations. The uncorrected signal firom a single 
channel appears in Figure 15 (top). Then calculated was lie average of 30 repetitions to 
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produce a teptesentation of the gradient activity. The averaging process removes noise 
uncorrelated with the imagiog gradients. Finally, this averaged signal was subtracted &pm 
the uncorrected signal, yielding the corrected signal appearing in the lower portion of 
Figure 15. After correction, the large amplitude gradient activity is removed completely 
5 and all that remains in the signal is the small, uncorrelated, noise fluctuation. 

The echo-planar imaging gradient artifact, which lasts for 23 ms and goes through 32 
sinusoidal oscillations in that period, is sampled with only four data points for these data. 
The digitized waveform is thus a very crude record of the actual gradient activity - well 
below the Nyquist frequency. In addition, the morphology of the gradient waveform is 
10 grossly different for each of the nineteen slices that are acquired in this TR period. This 
is due to the fact that the digitization is phase shifted s%htly with respect to the readout 
of each slice locatiotL There is also no amplifier "ring down" following the gradient 
pulses, as there is no higjh pass filter on the amplifier input. 

A similar experiment was performed, recording electrical potentials from the human 
15 scalp with a single differential electrode pair. Once the subject was placed into the 
imagtag system, the scanning protocol was similar to that used above: fourteen slice 
locations with a 3 second tr. As acquisition of each slice requires approximately 38 msec 
of gradient activity, the gradients can be expected to obscure about 20% of the EEG 
record As a consequence of several system non-idealities (such as the small inductance 
20 of the EEG leads, there will still be some electrical "ring-down." Figure 8 (bottom) 

shows 3 s of the EEG record during scanning, sampled at 200 Hz and triggered by the 
scaimer. Prominent gradient artifacts are present even after low pass filtering. Figure 16 
also shows a series of 20 successive 3 second traces after gradient artifact correction. 

EEG Energy Mapping 

• 2 5 The subject was scanned while he lay prone in the magnet with his eyes open. The raw 
EEG data were processed as outlined above^ first blanking the gradients and then 
removing the ballistocardiogram. The data were then submitted to a Fourier analysis and 
the energy at each of five pre-defined frequency bands was determined for each TR (that 
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is, for each image time point). Fot each frequency band, thetefote, it was possible to 
generate a separate time coutse. 

Shown in Figure 17A is the energy as a function of time in the alpha and theta bands for 
a four and a half minute period while the subject was at rest The energy levels in these 
5 bands are largely independent for the first minute and a half of this session and seem to 
co-vary for the latter portion of the recording. Figure 17B demonstrates the estimated 
time course of the BOLD signal response based on the (untested) assumption that the 
hemodynamic response related to the EEG is similar to that seen in activation smdies 
(e.g., (Cohen 1997; Cohen and DuBois 1999)). 

10 Using the latter as a reference function, then caloilated were correlation maps for each 
pixel location with the energy intensity in each previously defined firequency band using 
scanSTAT. The images were first subjected to a simple 9 pixel smoothing to reduce 
pixel noise. The resulting brain maps are shown in Figures 18A-C, with a separate map 
for the pixel correlation with each of the frequency bands. The images are remarkable in 

15 showing very high correlations, litde artifactual activation in white matter or CSF, and 
substantial symmetry, all of which suggest that the artifact content is low. Note also that 
the signal intensity in extra-sttiate areas shows a substantial signal decrease with increases 
in alpha level. 

Conchisions 

2 0 While diere have been numerous obstacles to the fusion of functional MRI and 

electroencephalographic data, the purely technical challenges of eliminating the cross 
contamination of the recorded signals seem to have yielded to rather straightforward 
engineering solutions of somewhat surprising simplicity. The tools outlined here are all 
both easy to construct and made fi:om inexpensive components. The complete sofibware 

2 5 solution can be made to run as a standalone unit using the LabView programming 

enviroimient The availability of high quality integrated instrumentation amplifiers such 
as the INA114 £tom Texas Instruments Corporation (Figure 14) makes the analog 
engineering uncomplicated as well 
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Certain factors affect the extent to which the data can be interpreted accurately. For 
example, the energetics of the spectral EEG signal are themsehres largely iinknown. 
Some components, such as the theta rhythm, have long been suspected to be driven by 
subcortical generators and might plausibly be associated with increased thalamic activity 
5 reflected as an increase in blood flow. The alpha rhythm, however, may well be intrinsic 
to the cerebral cortex (though influenced by thalamus) and is often interpreted as a sort 
of cortical resting state. What is clear, however, is that by continuing to study the 
relationship between scalp (and presumably braia) electrical potentials and magnetic 
resonance signal intensity (presumably indicating increases or decreases in cortical blood 
10 flow) it should be possible to understand better the physiological basis of the 
electroencephalogram (Schomer, Bonmassar etoL 2000). 

The signal processing methods developed to solve the problems in simultaneous EEG 
and £NfRI have a broad range of additidnal applications. For example, these methods can 
be used, without modification, to study electrical evoked potentials (Bonmassar, Anami et 

15 aL 1999) and their localization or, conversely, to yield a better temporal of the ''single- 
trial'' evoked MR response (Budkner, Bandettibi etaL 1996). The rliniral applications of 
fully integrated fMRI and EEG are substantial, the most immediate harvest will likely be 
m seizure source localization as an adjunct to oflier diagnostics used in surreal planning 
(Engel 2000). Interictal spikes are a common finding in the epileptic brain yet even when 

2 0 recorded fi:om the cortical sur&ce do not provide definitive and reliable source 

localization. Recendy, however, fast ripples have been reported in rat models and in the 
epileptic human brain, that are associated dosely with seizure foci (Bragin, Wilson et oL 
2000). A still speculative interpretation is that such activity represents rapid spontaneous 
action potentials that are ultimately propagated and observable as interictal spikes. If so, 

25 this gives added hope to the idea that localization by fMRI of interictal spikes might 

become a reliable means to plan resective brain surgery, as has been proposed already by 
several others (e.g., (Warach, Levin et aL 1994; Warach, Ives et oL 1996; Allen, Polizzi et 
aL 1998; Ramabhadran, Frost etaL 1999; Hoffmann, Jager etal. 2000)). 
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Aldiough EBG tecotded duxing MR imaging is patticulatly noisy, most of the atti&cts 
contended with here axe present in conventional EEG, though at lower amplitude. For 
example, the activity of the heart, particularly the electrocardiogram, is a contaminant 
that can be corrected in the same manner as the minimrzation of die ballistocardiogram. 
5 AC line noise is often present as well, and could be eliminated by using a tri^ered 
sampling approach, timing each block of samples to the power line oscillations — daese 
could be detected readily widi a phase-locked loop. 

Digital sampling based on precise synchroni2adon with well characterized noise sources 
(i.e., tri^ering on the basis of scanner gradient timing has applications that extend well 
10 beyond Magnetic Resonance Imaging. It is dear, for example, that a comparable 
approach could be used to remove noise sources associated with AC power line 
oscilladons; these include die removal of "hum" ficom digitally processed audio, light 
flicker £com fluorescent illumination, and other contaminants of digitally sampled signals. 

Example 3: Increasing Dynamic Range Available in Digitized Electro- 
15 encephalogtaphic Signals 

The scalp electdcal potentials used in the EEG contain both time varying and static (DC) 
components. Often the DC offset is much larger than the EEG, but it is seldom of 
^interest for clinical diagnostic purposes, as it contains essentially no information. 
However, it does cause tcouble as it increases the dynamic range needed to digitize the 

2 0 EEG signal For example, the EEG may be only a few microVolts ((xV), while potentials 
of several tens of milliVolts (mV) may exist betotreen electrodes, or as a result of the 
chemical electrode potential developed when the dectcode is in contact with the scalp. 
The depth of the signal digitization will be reduced by the ratio of the EEG to the DC 
potential: assuming a DC offset of 10 mV, an EEG signal of 10 jxV, and a 12 bit analog 

25 to digital converter (ADC), the 4096 different levels representable by the ADC will be 
reduced to only 4 levels for the EEG. Clearly this is unacceptable, as the quantization 
noise will dominate the recordings. 
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Fot this reason, conventional EEG am^lifieis ate equipped with AC-coupled Qiig^ pass) 
inputs, usually a capacitor separating die output of die first stage amplifiers fi:om the 
input to die ADC. As the firequencies of interest in EEG can be quite low» die inputs will 
typically have time constants of several seconds, allowing fluctuations of 1 H2 or so to 
5 pass widiout significant attenuation. One consequence of diis AC coupling is diat it 

creates a time constant for signal recovery if die input saturates and, because diese filters 
must pass very low firequencies, die setding time for die analog signal to come back to 
die center of its nominal range can be quite long. For conventional EEG recordings diis 
is acceptable, as laige DC shifts are not fi:equent 

In £MRI, die recovery time associated with AC coupling is a partictdar problem, as the 
gradient-induced artifacts can be large enougjh to bring die input stages into saturation. 
When diis occurs, die amplifiers diat follow can be pinned at either die positive or 
negative supply rails for several milliseconds. When die gradients cease, die amplifiers 
may require significant setding time such diat die gradient artifact substantially outlasts 
die gradient event 

A recent report for example, firom Lovblad, et al., reported that even when using a 
reduced gradient activity BURST (Hennig and Hodapp 1993) sequence, "The EEG 
could not be interpreted dviring die artifact caused by the excitation pulse, but the 
recording becomes readable in less dian 1 second (approximately 100 msec) after 
2 0 completion of die BURST" (Lovblad, Thomas et al. 1999). Thus, if they were to use, for 
example, an eight slice acquisition with a 2 second TR (repetition time), die EEG signal 
would be obscured more dian half of die time by eidier gradient activity or amplifier 
recovery. 

To remove the DC offset, and to reject saturation problems, the invention provides a 
2 5 circuit tiiat behaves as illustrated in Figure 19. 

In this circuit, the high firequency common mode artifacts typically present in EEG 
recorded in die MRI environment are first attenuated using passive components before 
being differentially amplified. The differential amplification furdier reduces common 
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mode artifacts. The data are then presented to a hig^ pass ('^anti-aliasing") filter to both 
ajttenuate the arti&cts and to minimize the possibility of aliasing in analog to digital 
conversion. The final stage of this circuit performs the dtial role of amplifying the signal 
further and removing any DC bias. This is done by comparing the signal differentially to 
5 a DC reference value sampled fi:om the signal itself The sampled signal is stored under a 
command here shown diagrammatically as a switch connected to ground. 

Figure 14 shows one preferred implementation, in which the first hig^ pass filter consists 
of a resistor-capadtor (RC) network 20. The fibcst differential amplifier is an integrated 
10 instrumentation amplifier 28. The active high pass filter is a 24 dB Butterworth filter. The 
sample/hold circuit 24 stores the DC signal as chaige across a low-leakage capacitor 37 
with a compensation circuit for the non-zero input offset current found in the amplifier 
that drives it The final differential stage is a standard operational amplifier 34. 

With more specific reference to Figure 14, a pair of electrodes is connected to the 
differential input amplifier stage 20, bypassed at high fiequencies by parallel capacitance. 
The input amplifier 28 itself may be an integrated circuit, such as an INAl 14 firom Burr- 
Brown corporation. A pair of matched 5 megohm resistors supply bias current to the 
integrated circuit. The output firom the amplifier is used as an input to a sample and hold 
stage 24, which utilizes an integrated circuit 36 such as die LF398 from National 
Semiconductor corporation, which stores the DC voltage across a capacitor 37. A series 
resistance of 200 kQ provides a time constant of 5 seconds to the sample and hold stage. 
The output of the sample and hold stage is appHed to the offset reference pin of the 
differential amplifier component. An offset trim adjustment 40 is provided to adjust for. 
small static DC offsets in the integrated circuit 36. An input 39 is provided to sample and 
hold device 36. 

The DC corrected output is presented to the iuput of a low pass filter 22, that utilizes 
operational amplifier integrated circuits 30 and 32, such as the TL072 firom Texas 
3 0 Instruments corporation. With the component values indicated in Figure 14, the low- 
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pass filter is implemented in a Chebyshev configunation to ptovide. apptoximately 30 
decibels/ octave of attenuation at fi:equendes above the selected pass band. A final gain 
stage 26 includes a gain trimTning adjustment 38 to match the gains actoss channels. 

5 Example 4: Futthet Reduction of Sensitivity to Electrical Interference From 
Ongoing Scanning Activity 

In an additional preferted embodiment, diagrammed in Figure 20, the differential 
recording input stages are supplemented by a shield driver circuit that furdier reduces the 
sensitivity to electrical interference firom the ongoing scanning activity, or other common 

10 mode signals in the leads. In this embodiment, the output firom die initial fltnplififtrR is 
coupled dirough an isolation amplifier to reduce the possibility of circulating currents 
firom the RF pulses saturatir^ later components of the amplifiers (Figure 20). The DC 
offset signal is first detected digitally, then held in a digital latch, and returned via a digital 
to analog convertor as an offset correction. This offset correction is also coupled 

15 throii^h an isolation amplifier. In tiiis embodiment, the analog output is coupled optically 
tbrbugh the RF shielded room to minimise any problems of corruption of the MR signal 
through Actors such as ground loops. Provision is tnade for the recording of multiple 
channels of EEG data through the use of multiplexing drcuitty. 
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The foregoing description of preferred embodiments of the invention has been 
presented for the purposes of illustration and descrqition. It is not intended to be 
exhaustive or to limit the invention to a precise form disclosed. Many modifications and 
5 variations are possible in light of the above teaching. It is intended tiiat tiie scope of the 
invention be limited not by this detailed description, but rather by the rW^m^ appended 
hereto. 



58 



02/13689 



PCTAJSOl/25480 



WHAT IS CLAIMED IS: 

1 . A method of tedadng contanmiation of electncal signals tecotded in the 
presence of repeated intetfetence contamination comprising: 

(a) obtaining an dectdcal signal^ wherein the electdcal signal was recorded in 
die presence of a contaminating signal; 

(b) detecting a timing signal that occurs at a fixed time point during the 
electrical signal relative to the onset of the contaminating signal; 

(c) digiti2ing the electrical signal, wherein the digitizing begins with the 
timing signal; 

(d^ analyzing a plurality of digitized electrical signals, wherein the electrical 
signals are synchronized with respect to the timing signal, to obtain an 
estimated contaminating signal; and 

(e) subtracting tiie estimated contaminating signal from the digitized 

electrical signal, thereby reducing contamination of the electrical signal 

2. The method of claim 1, wherein the analyzing to obtain the estimate of the 
contaminating signal comprises averaging the electrical signals. 

3. The method of claim 1, wherein the analyzing to obtain the estimate of die 
contaminating signal comprises calculating a weighted average of die electrical 
signals. 

4. The method of claim 3, wherein the estimate of the contaminating signal is 
biased towards recent events. 

5. The method of claim 4, wherein the estimate of the contaminating signal 
comprises a first sum consisting of the ndi electrical signal added to a scalar 
multiple, w, of a prior estimate of the contatninating signal, divided by a second 
sum obtained by adding a series 1+w^w^+w^ ... + ... w^. 
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6. The method of clattn 1, 2 ox 3, \dieteui the electdcal tecording comptises an 
decttophysiological signal. 

7. The method of claim 1 , wherein die estimate of the contaminating signal is 
multiplied by a scalar prior to the subtracting of step (e). 

5 8. The mediod of claim 4, wherein die electrophysiological signal comprises an 
electroencephalographic recording, an electromyelographic recording an 
decttocatdiographic recording or a measure of galvanic skin resistance. 

9. The mediod of any one of claims 1-8, wherein die interference comprises 
interference arising fi:om inductively coupled magnetic fidds. 

10 10. The method of any one of claims 1-9, wherein die interference comprises 
interference arising £t:om alternating current (AC) line noise. 

11. The mediod any one of claims 1-10, wherein die digitizing is performed at a 
sampling rate bdow die Nyquist rate for die contaminating signal 

12. The mediod of any one of claims 1-1 1, wherein die dectrical signal obtained in 
15 step (a) is passed dirough a low pass filter prior to the digitizing, at a frequency of 

approximatdy one half of die frequency at which die electrical signal is sampled. 

13. The method of claim 12, wherein the low pass filter passes signal fisquendes of 
less dian about 200 Hz. 

14. The method of any one of claims 1-13, wherein the method is performed 
2 0 concurrentiy widi Magnetic Resonance Imaging of die subject. 

1 5. The mediod of claim 14, wherein die dectrical signal comprises an 
dectrophysiological signal and the contaminating signal comprises gradient 
activity. 
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16. The method of claim 15, wherein the electrical signal comprises an 
electrophysiological signal and the contaminating signal comprises radio 
firequency transmitter activity. 

17. The method of any one of rlaims 1-1 6^ wherein the digitizing is performed at a 
5 rate of about 200 to about 5000 samples per second. 

1 8. The method of claim 1 , further comprising removing a DC offiset ficom the 
electrical signal by analog subtraction prior to the digitizing. 

19. The method of claim 18, wherein the DC offset is measured and subtracted from 
the electrical signal using a differential amplifier. 

10 20. The method of claim 19, wherein the DC offset is measured by analog to digital 
conversion, and averaged over a time period long compared to lowest 
frequencies of interest in the electrical signal 

21 . The method of claim 20, wherein the time period is approximately 10 times 
longer than the lowest frequencies of interest in the electrical signal 

15 22. The method of claim 20 or 21, wherein the analog subtraction comprises 

converting the averaged signal to an analog voltage and electrically subtracting 
the averaged signal from the electrical signal through differential amplification. 

23. The method of claim 19, wherein the DC offset is measured in an analog 
integrator having a time constant long compared with lowest frequencies of 

20 interest in the signal 

24. The method of claim 23, wherein the time constant is approximately 10 time the 
lowest frequencies of interest. 

25. The metiiod of daim 7, wherein die electrophysiological recording comprises an 
electroencephalogram and is recorded concurrendy with magnetic resonance 

25 image acquisition. 
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26. The metbod of claim 25, herein the electtophysiolo^cal tecordings ate used to 
infotm interpretatioiis of the magaetic resonance images. 

27. The methods of claim 25 ot 26, wherein the electrophysiological recordings are 
used in a statistical analysis of change in intensity of the magnetic resonance 

5 signal. 

28. The method of any one of rlaitna 25-27, further comprising determining a 
correlation between change in intensity of the nmgnetic resonance signal and a 
feature of the electrophysiological recording. 

29. The method of claim 28, wherein the correlation is used to make statistical 
10 images that represent an association between the electrical signals and the 

intensity of the magnetic resonance signal intensity. 

30. The mediod of claim 28, wherein die feature of the electrophysiological 
recording comprises a dme course of signal intensity change in defined frequency 
bands contained in the electcophysiological recording. 

15 31. The mediod of claim 30, wherein the defined firequency bands correspond to 
standard ranges used for clinical interpretations of the electroencephalogram. 

32. The method of daim 31, wherein the standard ranges are selected firom the group 
consisting of ficom 0 to approximately 4 Hz (the Delta band), from approximately 
4 to approximately 8 Hz (the Theta ban(^, ficom approximately 8 to 

20 approximatdy 12 Hz (the Alpha band), fi:om approximately 12 to approximately 

30 Hz (the Beta banc^, and firom approximately 30 Hz and greater (the Garmna 
band). 

33. The method of any one of claims 25-32, further comprising convolving the time 
course of the electrophysiological signal vnih an estimate of the magnetic 

2 5 resonance hemodynamic impulse response fimction. 
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34. A medaod of reducing magaetic intetfexence diiidng 

tecotding £com a subject by measumig an electdcal potential difference between 
a pair of electrodes, wherein the pair of electrodes conununicate with a 
differential amplifier via electrical connections, the method comprisiDg twisting 
5 the electrical connections together, thereby reducing magnetic interference. 

35. A method of reducing magnetic interference during electrophysiological 
recording £com a subject by measuring an electdcal potential difference between 
a pair of adjacent electrodes, wherein each electrode comprises two leads, the 
method comprising twisting each lead together with a lead of an adjacent 

1 0 electrode, thereby reducing magnetic interference. 

36. The method of ckim 34 or 35, wherein the decttophysiological recording 
comprises an electroencephalographic recording. 

37. The method of any one of claims 34-36, wherein the method is performed 
concurrendy with Magnedc Resonance Imaging of the subject 

15 38. An apparatus for processing digiti2ed electrical signals in the presence of a 
repeated contaminating signal comprising: 

(a) a signal processor 90 adapted to receive a recording of an electrical signal; 

(b) a detector 92 adapted to detect a timing signal that occurs at a fixed time 
point dudng an dectrical signal relative to the onset of a contaminating 

20 



(c) a signal accumulator 94; and 

(d) a processor 96 adapted to subtract averaged waveforms from an dectrical 
signal. 
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